Detecting an exciton crystal by statistical means 
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We investigate an ensemble of excitons in a coupled quantum well excited via an applied laser 
field. Using an effective disordered quantum Ising model, we perform a numerical simulation of 
the experimental procedure and calculate the probability distribution function P(M) to create M 
excitons as well as their correlation function. It shows clear evidence of the existence of two phases 
corresponding to a liquid and a crystal phase. We demonstrate that not only the correlation function 
but also the distribution P(M) is very well suited to monitor this transition. 

PACS numbers: 73.21. La, 72.70.+m, 73.23.-b 



The exciton is a very fascinating composite particle, 
which can be generated and investigated in specially de- 
signed semiconductor heterostructures - the bilayer sys- 
tems. In its simplest incarnation it is a bound state of 
an electron and a hole and is thus of bosonic nature. If 
the system size of such a compound is small it is natu- 
ral to ask, whether a Bose-Einstein condensation (EEC) 
of such structures is possible [1]. On the other hand, 
also the possibility of a Cooper-pair-like ground state [the 
conventional Bardeen-Cooper-Schrieffer (BCS) supercon- 
ductivity] has been considered in a number of works [2] . 
Recent experimental progress in the field of electronic 
bilayer systems allowed for a study of all these different 
fascinating possibilities O H] . 

In most experimental realizations the holes and elec- 
trons are spatially separated so that every such indirect 
exciton has a relatively large dipole moment. It turns 
out that at some intermediate density, at which the BEC 
condensation is prohibited, the dipole interactions let the 
excitons see each other. If the correlations are strong 
enough even a long-ranged ordering of Wigner crystal 
type is possible [5^, 'BJ . A detection of such kind of crys- 
talline structure is very difficult though. The excitons 
themselves are usually generated with the help of laser 
fields and they are rather fragile with respect to irradia- 
tion. Thus traditional spectroscopic techniques are very 
difficult to apply and one needs alternative methods [5]. 
One such approach is based on the knowledge of the ffist 
order correlation function, the measurement of which was 
very recently reported in Ref. [4]. 

Here we propose an alternative statistical method of 
detecting and analyzing the properties of the exciton 
crystallization phenomenon and discuss its predictive 
power. A typical experimental cycle would start with 
the generation of excitons via a laser pulse in a coupled 
quantum well structure (e.g. GaAs/AlGaAs heterostruc- 
ture) . After that the number of excitons M is measured 
by their recombination. Conducting a large number of 
cycles one gathers the statistics of M [7l[8]. The fact that 
for a given M a regular arrangement of the excitons on a 



lattice minimizes their interaction energy should be vis- 
ible in the probability distribution P{M). Although the 
correlation function possesses a higher predictive power, 
we shall show below that the crystallization can even be 
seen in P{M)^ which is accessible by much less effort. 
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FIG. 1. (a) shows the typical experimental situation: excitons 
are created in a coupled quantum well via an applied laser field 
that excites electrons from the valence band into the conduc- 
tion band and thus forms excitons. Possible exciton states 
in a typical coupled quantum well such as GaAs/AlAs are 
shown in (b) . These coupled quantum wells have two conduc- 
tion bands (F and X) originating from different points in the 
Brillouin zone and allow for both the formation of direct and 
indirect excitons, however, with different binding energies. 

We assume that before every measurement cycle the 
system consists of TV valence band electrons, located at 
random positions r^. We model each electron by a two- 
level system consisting of an unexcited state and an exci- 
ton state. A laser field with frequency ujl excites the sys- 
tem. It is detuned from the actual transition frequency 
between these states by A. By means of the laser inten- 
sity the Rabi frequency cj, which accounts for transitions 
from ground to exciton state, can be changed. In a typ- 
ical experiment [9 one uses short laser pulses with high 
intensity, for which the Rabi frequency is of the same 
order as the exciton binding energy [10] . A circularly po- 
larized light beam (polarization a) at a suitably chosen 
frequency can create direct and indirect neutral excitons 
with well defined spins a and a = — cr by exciting an 
electron from the valence band of the same or the other 
semiconductor in the coupled quantum well (layer), re- 
spectively (see Fig. fl]) [S^illj. 



In the following we will focus only on the indirect ex- 
citons, the particle and hole of which are located in dif- 
ferent layers. The formation of direct excitons can be 
neglected, since either their formation can be suppressed 
by an appropriate choice of the excitation frequency cjl, 
or one can simply wait long enough - their recombina- 
tion time is much shorter. The indirect excitons possess 
dipole moments d perpendicular to the layers d = eD^ 
where D is the interlayer separation, see Fig. Ilia). 

We treat the exciton as a two level system [12] with 
electrons either in the valence band or excited to an ex- 
citon [13]. We assume the levels to be sharp neglecting 
effects from the Fermi distribution of the separate bands. 
We can do so in the limit of a large detuning of the laser 
from the resonance [10]. The electrons interact with the 
laser light and with each other through the dipole interac- 
tion when they form excitons. The velocity distribution 
of excitons in coupled quantum wells [9] can be tuned by 
efficient cooling [4 and the application of a perpendicu- 
lar magnetic field, which gives rise to a higher effective 
mass [HI [14] . With both the laser pulse duration [15] and 
the experimental measurement taking r < 0.1 ns and the 
resulting velocities vex ^ 10 m/s the displacement of a 
single exciton is vex'T' < 1 nm. The typical separations 
of the excitons one usually encounters are of the order 
100 nm so that we can assume the excitons to be fixed in 
space [16] . The size of an exciton is estimated via its Bohr 
radius as ~ 20nm, see e.g. [17]. Since the electrostatic 
properties do not depend on the details of exciton states 
we model the system as a randomly arranged interact- 
ing ensemble of spin 1/2 sub-systems each representing 
a single electron/exciton. Hence, the Hamiltonian reads 
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where ax,z denote the Pauli matrices. This equation 
uses the rotating wave approximation when describing 
the light-matter interaction as in Ref. [19]. It neglects 
all terms oscillating with frequencies ujl and higher. 
From here on we measure all energies in units of A en- 
tailing u = cj*/A. The interaction strength is there- 
fore measured in units of AL^, leading to the defini- 
tion C = C*/AL^. The ground state particle density 
n = N/V ^ where V is either a one-dimensional interval 
or a two-dimensional square, is given for all plots. 

The model ([l]) can be interpreted as a (generalised) 
spin-1/2 Ising model. The Rabi frequency uo* and the 
detuning A correspond to magnetic fields in x- and z- 
direction. The third parameter C* indicates the strength 
of the effective interaction between the excitons. We 
should note that a similar model has originally been ap- 



plied to interactions between Rydberg atoms, where a 
similar crystal-like phase exists [HI [20] . In the case of 
Rydberg atoms, induced dipole moments gave rise to 
van-der-Waals interactions. In the case of excitons, in 
contrast, dipole-dipole interactions between the excitons 
lead to a stronger dependence of the distance, ex \r\~^ . 

We are only interested in the case A > 0, because 
otherwise it is energetically not favorable to produce ex- 
citons. We typically use a large detuning from the tran- 
sition frequency in accordance with experimental studies 
[9 , so that A is larger but still comparable in magni- 
tude to the Rabi frequency. A typical laser field has an 
excitation frequency of about 1 eV. The excitons have a 
dipole moment oriented perpendicular to the plane. In 
this case C* = e^D'^ /e. For the dielectric spacer between 
the top and bottom layer we assume e = 12.9eo being a 
typical value for GaAs and D = 11.5nm as put forward 
by [21 . In the numerical simulations we take the length 
L of the simulated square in 2D to be ~ 200 — 500nm. 
For such and larger system dimensions we did not detect 
any sizeable finite size effects. 

The numerical procedure emulates the experimental 
process by initially generating a random distribution of 
N electrons [spins in Eq. (fl])] in a given one- or two- 
dimensional volume with open or periodic boundary con- 
ditions. Then the corresponding Hamiltonian matrix is 
set up. For large N the size of this matrix is reduced 
by the truncation of the Hilbert space, which is done by 
taking into account only k basis states with the smallest 
diagonal elements in the Hamiltonian matrix. We sys- 
tematically checked that all the results do not depend 
on k. In the next step the eigenvector corresponding to 
the smallest eigenvalue, the ground state (GS), is calcu- 
lated by matrix diagonalization [22 . It is given as a lin- 
ear combination of the previously mentioned basis states 
and therefore enables an efficient evaluation of different 
observables. Besides the number of excitons (which can 
be non-integer since the ground state is generally a su- 
perposition) the pair correlation function ^(r), which is 
closely related to the density distribution of excitons, can 
be computed in the following way: we divide the possi- 
ble range for distances between electrons in our system 
into equidistant bins and measure the distance between 
each pair of electrons to assign it to a certain bin. For 
each pair the squares of coefficients from the representa- 
tion of the GS as a linear combination are summed over 
those states in which the particular pair of electrons is 
excited. The sum over multiple random arrangements 
then produces the correlation function. 

We start by considering a ID system with different in- 
teraction strengths. Already for weak interactions we see 
that the correlation function is zero for r < Rb, where 
Rb can be interpreted as the blockade radius, see Fig. [2J 
This effect is due to strong dipole field in vicinity of an ex- 
citon which suppresses the excitation of additional ones. 
As expected Rb grows for increasing interaction strength. 
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FIG. 2. 
and nL 



Correlation function for the ID case for cj = 0.1 
= 25 for different interactions strengths. The open 



circles correspond to C = 7 • 10 , the black triangles refer to 
C — 7 ■ 10~^ and the black squares are for C — 7 ■ 10~^ . 
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FIG. 3. Variance of the mean number of excitons for uo — 0.1, 
and nL — 25. The straight lines and the arrow illustrate the 
phase transition. The inset shows the excited fraction of exci- 
tons, which is bounded by 1 in the limit of weak interaction. 



Simultaneously an emergence of peaks of higher order is 
observed. For instance, at C > 10""^ the second order 
maximum is already very pronounced and the pair cor- 
relation function resembles that of a liquid. Finally, at 
C > Ccrit ^ 10~^ g{r) drops to almost zero between the 
peaks, which we attribute to the crystalline ordering of 
the excitons which is caused by 'bubbles' around each 
exciton in which no further excitation is possible. 

It is instructive to analyze how this evolution of the 
system is reflected in the distribution P{M). The sim- 
plest distribution parameters are the mean (M) and the 
variance a^ = (M^) - (M)^, which we plot in Fig.js] Es- 
pecially the latter quantity significantly changes its be- 
haviour in vicinity of the critical interaction strength. 
For weaker C there is only very little change in a^ while 
for C > Ccrit it is described by a power law. This dras- 
tic change is due to formation of the dipole crystalline 
order. Thus, just by measuring P{M) one can identify 
the interaction strength at which the crystallization takes 
place. 

The system under consideration features long-ranged 
interactions so that the coordination number for every 



particle is effectively infinite. That is why it is reason- 
able to expect that the actual dimensionality of the sys- 
tem should not play any significant role in the physical 
picture. In order to substantiate this argument we have 
extended our analysis to 2D systems. Here one has to 
work with a larger number of particles. As a result one 
is confronted with considerably longer computation times 
(they grow exponentially with the number of particles) . 
Nonetheless one is able to perform the same numerical 
program as in a ID case. 

In Fig. [4] we plot the correlation function in 2D case, 
which is generated for a system with periodic boundary 
condition and 30 spins. For weak interactions we again 
obtain a 'flat' correlation function, which is essentially 



uniform for r > Rb- For C > C^ 



10- 



(r) starts 



to develop additional maxima, which is a signature of 
the crystallization onset. As in the ID case we would 
like to simultaneously monitor the evolution of P{M)^ 
see Fig. [5] for the plot of the average and variance. The 
variance a^ [as well as g{r)] shows a clear transition at 
C ~ 10~^. The change of the power law decay of M 
is not as abrupt as in another possible realisation of a 
Wigner crystal though ^^. The remarkable change of 
the variance as a function of C shows that also simple 
statistical analysis of the exciton number allows to iden- 
tify the phase transition to an exciton crystal. 

In conclusion we have investigated the formation and 
detection of an exciton crystal in a typical semiconduc- 
tor environment using exact numerical diagonalization 
and approximative descendants of this method for one- 
and two-dimensional systems. Both the pair correlation 
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FIG. 4. Pair correlation function for uj = 0.1, C = 8 • 10~^ 
at nL^ = 25 (black squares). The blockade radius as well 
as the first two maxima of the curve are clearly visible. The 
double peak of the second maximum really is a single peak 
and only splits due to commensurability with the periodic 
boundary conditions in the sample shown. The data points 
for the largest distances are less reliable than the rest since 
only very few runs contribute to these points. We compare 
our result to the case C = 2 • 10~^ where no additional peaks 
and only the blockade radius is present (open circles). The 
inset shows the histogram of the number of excitons with a 
Gaussian fit for the same parameters as for g{r). 
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FIG. 5. Excited fraction (left panel) and variance (right panel) of a 2D system (this is a direct analogon of Fig. 3 of the ID 
case). Different data sets correspond to different particle numbers at equal density. The coding is as follows: A^ = 10 (yellow, 
solid squares), N = 15 (orange, solid triangles). A/" = 20 (red, open circles), N — 25 (purple, open squares) and A/" = 30 (blue, 
sand clocks). One observes a very good convergence for increasing N . Taking the edge length of the system with A^ = 30 as a 
reference one can map all the curves on it by rescaling of L, which is done in the insets by streching the abscissa by a factor of 
{Ln=3o/Ln) ' 



function and the simpler statistical data show signatures 
of this phase transition. The blockade phenomenon dis- 
tinguishes the crystal- and liquid-like behavior from the 
BEG phase, whereas additional peaks in g{r) show the 
phase transition from the liquid to the crystal phase. 
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